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Abstract. The purpose of the paper is three-fold: (a) we prove that every sequence which is a multidimen- 
sional sum of a balanced hypergeometric term has an asymptotic expansion of Gevrey type-1 with rational 
exponents, (b) we construct a class of G-functions that come from enumerative combinatorics, and (c) we 
give a counterexample to a question of Zeilberger that asks whether holonomic sequences can be written as 
multisums of balanced hypergeometric terms. The proofs utilize the notion of a G-function, introduced by 
Siegel, and its analytic/arithmetic properties shown recently by Andre. 



1.1. Balanced multisum sequences are of Nilsson type. The purpose of the paper is three- fold: 

(a) we prove that every sequence which is a multidimensional sum of a balanced hypergeometric term 
has an asymptotic expansion of Gevrey type-1 with rational exponents, 

(b) we construct a class of G-functions that come from enumerative combinatorics, and 

(c) we give a counterexample to a question of Zeilberger that asks whether holonomic sequences can be 
written as multisums of balanced hypergeometric terms. 

The proofs utilize the notion of a G-function, introduced by Siegel, and its analytic/ arithmetic properties 
shown recently by Andre. Let us begin by introducing the notion of a (balanced) multisum sequence. 
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Definition 1.1. A (balanced) multisum sequence (a n ) is a sequence of complex numbers of the form 

(1) a n = 2J *n,fc 

fc6supp(t ni .) 

where t is a (balanced) term and the sum is over a finite set that depends on t. 

Definition 1.2. A term ^ in variables (n, k) where k = . . . , k r ) is an expression of the form: 

r J 

(2) t^c^n^n^'^ 

»=i i=i 

where Cj € Q for i = 0, . . . , r, ej — ±1 for j — 1, . . . , J, and A,- are integral linear forms in the variables 
(n, fc) such that for every n £ N, the set 

(3) supp(t n ,.) := {k£l7\ Ajin, k) > 0, j = 1, . . . , J} 

is finite. We will call a term balanced if in addition it satisfies the balance condition: 

J 

(4) E e ^ = °- 

3=1 

For example, the Apery sequence (see |vdP] ) is a balanced multisum given by: 




Multisum sequences appear frequently in enumeration questions; for numerous examples, see |St2[ IFS] . A 
key problem is to study the asymptotics of a (balanced) multisum sequence. This is a classical problem that 
has been discussed by several authors for decades, see |BTllFSUSt2[|WP] . Parsing through the literature, in 
numerous examples of balanced multisum examples, a certain rationality of the leading exponents of n was 
found by accident, with no explanation. Understanding this rationality lead to the results of our paper. 
To explain this rationality, let us introduce one more definition. 

Definition 1.3. We say that a sequence (a„) is of Nilsson type if it has an asymptotic expansion of the form 

(6) J2 ^ n n a (log(n) ff x ^ 

X,a,0 

where the summation is over a finite set of triples (A, a,/3), A £ Q is an algebraic number, a £ Q, f3 € N, and 
f\,a,p(z) is Gevrey-1 power series, i.e., the coefficient of z k in f\. a ,p(z) is bounded by C n for some C > 0. 

Now we can state our first result which does not seem to be covered by the existing literature on asymptotic 
expansions of sequences. 

Theorem 1. Every balanced multisum sequence is of Nilsson type. 

1.2. G-functions. The proof of Theorem [T] utilizes the notion of a G-function, introduced by Siegel in 
[Si] with motivation being arithmetic problems in elliptic integrals, and transcendence problems in number 
theory. For further information about G-functions and their properties, see |An| IBo| IDGS[ I5 H IToj . 

Definition 1.4. We say that series G(z) = J2^Lo a " 2 " i s a G-function if 

(a) the coefficients a n are algebraic numbers, 

(b) there exists a constant G > so that for every n € N the absolute value of every conjugate of a n is 
less than or equal to G™, 

(c) the common denominator of ao, • ■ ■ , a n is less than or equal to G™, 

(d) G(z) is holonomic, i.e., it satisfies a linear differential equation with coefficients polynomials in z. 
Remark 1.5. In [Anj . Andre calls a series that satisfies (a)-(c), an arithmetic Gevrey-0 series. 
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Our next theorem is a construction of G-functions from balanced terms. Notice that it is easy to generate 
examples of balanced terms; see for example the Apery sequence above. 

Theorem 2. If (a n ) is a balanced multisum sequence, the generating series G(z) = X)^Lo a « z ™ * s a G- 
function. 

G-functions can be constructed by arithmetic, or by geometry. Theorem [2] offers a new construction. For 
further discussion, see Sectional 

The proof of Theorem [2] also gives the following result, which may be of interest to enumerative combi- 
natorial problems that often lead to sequences of integers. 

Theorem 3. The generating series of an integer-valued, exponentially bounded holonomic sequence is a 
G-function. 

1.3. Holonomic sequences are not balanced multisums. We now come to the third part of the paper, 
which compares multisum and holonomic sequences. 

Definition 1.6. A sequence (a n ) is holonomic (i.e., D-finite in the sense of [St lj ) if it satisfies a linear 
recursion relation with polynomial coefficients. In other words, there exist d € N and Pj(n) G Q[n] (where 
Q denotes the set of algebraic numbers) for j = 1, . . . , d, so that for every n £ N we have: 



The following is a fundamental theorem of Wilf-Zeilberger. 

Theorem 4. |WZ1 IZe] Every multisum sequence is holonomic. 

The above theorem has a constructive proof with several computer implementations, see |PR| IPWZ] 
and [WZ . The converse was widely accepted as a reasonable conjecture, communicated to the author by 
Zeilberger. Our goal is to give a counterexample, and give an obstruction for the converse to hold. 

Theorem 5. Consider the holonomic sequence (a n ) defined by 



with initial conditions oq — 0, a± = 1. Then, (a n ) is not a balanced multisum. 

To understand why the converse to Theorem |4] fails, and why the example given by (JSj) is not pathological 
(but rather typical), let us look at the asymptotic expansion of an exponentially bounded holonomic sequence. 
It follows from Birkhoff-Trjitzinsky and Turrittin that a holonomic sequence is almost of Nilsson type, i.e., 
it satisfies ([6]) where the exponents a are algebraic, but not necessarily rational numbers. This also typically 
happens in the analysis of linear ODE, where the above exponents are known as Frobenius exponents. For 
more details, see Section 12.41 On the other hand, balanced multisum sequences have rational exponents 
according to Theorem [TJ This proves and explains Theorem [5] 

1.4. Acknowledgement. The author wishes to thank Y. Andre and D. Zeilberger for enlightening conver- 
sations, guidance and encouragement. 



2.1. Proof of Theorem [2j Let us begin with the following alternative presentation of a balanced term. 
Lemma 2.1. Every balanced term i can be written in the form: 



(7) 



P d (n)a n+d + . .. P (n)a n = 0. 



(8) 



(2n + l)a n+2 - (7n + ll)a n+ i + (2n + l)a n = 0. 



2. Proofs 



(9) 




where Ci G Q for i = 0, . . . , r, ej = ±1 for j = 1, . . . , J, and Bj, Dj are integral linear forms in the variables 
(n,k). 
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Proof. Consider a balanced term t given by ^ , where the linear forms Aj satisfy the balance condition (01 . 
Let J ± = {j e J \ej = ±1} and consider the linear form A(n, k) defined by: 

A(n,k) = Aj(n,k) = V" Aj(n,k), 
jeJ+ jeJ- 

where the second equality follows from the balance condition. Then, multiply and divide the balanced term 
by A(n, k)\, and rearrange the factors into a ratio of multibinomial coefficients as follows: 



Y[(Mn,ky.y> = 



Y\ ]GJ+ A 3 {n,k)\ A{n,k)\ 
U je j- Aj(n,k)\ A(n,k)\ 

I A(n,k) \ 

va 3 \jej-) 



A(n,k) \ 



\Aj \j£j+) 

Now, write the multibinomial coefficient as a product of binomial coefficients. The result follows. □ 

The next lemma from number theory is well-known (see |vdP|, p. 198] and also [Si]) and follows easily from 
Chebytchev's theorem. Below, 1cm denotes the least common multiple. 

Lemma 2.2. [^dPj [Si] Th ere exists C > so that 



(10) lcml (?],...,("] )<cr 



for all n G N. 

Proof. Let ord p m denote the maximal power of a prime number p that divides a natural number m. Then, 
for every natural number n and integer a and b with < b < a < n and every prime number p we have: 



ordp 







'logo" 


C) 


< 








.logp. 



ord p 6 < ord p lcm(l, . . . , a) — ord p 6 < ord p lcm(l, . . . , n). 



Thus, 



lcml ["),...,(")) <lcm(l,...,n) 



On the other hand, it is known that 

loglcm(l, . . . , n) = 0(n). 

For a detailed discussion, see |vdP[ p. 198]. □ 

We are now ready to give the proof of Theorem [21 

Proof, (of Theorem [5) Fix a balanced term t^/. as in and the corresponding sequence (a„) of (fT]). We 
will show that conditions (a),(b),(c) of Definition 11.41 are satisfied. Condition (a) is obvious. 

Using (?) < 2 a , and the fact that the set is a subset of [— Kn,Kn] r n Z r for some K > 0, Equation 
© implies that there exists a constant C > so that 

IWI <° n 

for all (n, k) and for all complex conjugates of t n ^. Summing up with respect to k in Equation ([T]), and 
using the fact that the summation set has polynomial size in n, it follows (after possibly enlarging C) that 

Knl <C" 

for all n > and for all complex conjugates of a n . This proves condition (b) of Definition 11.41 

Condition (c) follows from Equation ([9]), Lemma [2721 and the fact that the summation set ([3]) is bounded 
poly normally by n. 
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2.2. The local monodromy of a G-function. In this section we will make little distinction between a 
convergent power series, its analytic continuation, and the corresponding function. Recall that a power 
series is holonomic if it satisfies a linear differential equation PG(z) = where P G Q(z,d/dz) is a linear 
differential operator with coefficients in Q[z]. By the theory of differential equations (see for example flnllO] ). 
a holonomic function has analytic continuation as a multivalued analytic function in C \ A, where A is a 
finite set of algebraic numbers. The following Theorem follows from a combination of results of Katz, Andre 
and Chudnovsky; see |Ka[ IAn[ ICC] and also |C-L) for a detailed exposition. 

Theorem 6. [An, Ka. CC The local monodromy T of a G-function around a singularity is quasi-unipotent. 
In other words, 

(11) (T r -l) s = 
for some nonzero natural numbers r and s. 

It follows that the local expansion of G(z) around a singularity A € A is a finite sum of series of the form: 

(12) ]T c a ^{z - A)«(log(z - X)fh a ,p(z - A) 

where a € Q, [3 S N, c a .p S C and h at p{w) are convergent germs at w = 0. In fact, Andre shows that 
ha,p{w) are G- functions themselves. Said differently, the G-functions that come from arithmetic 

(a) are regular holonomic (i.e., the power series h a ,p{w) above are convergent at w — 0), and 

(b) have rational exponents (denoted by {a} above). 

On the other hand, the generating series of a generic exponentially bounded holonomic sequence (a„) will 
not be in general regular holonomic, nor will they have rational exponents. This explains Theorem [5] 

Remark 2.3. Power series of the form (I12p are known in the literature as Nilsson series; see [Nij . 

2.3. The Taylor series of a G-function and Theorem Q3 The following lemma is a well-known appli- 
cation of Cauchy's theorem; see for example [Ju| Thm.A]. 

Lemma 2.4. [HI Thm.A] If a E C\N, (3 e N, and 

(l-zr(log(l-z)y 3 = J2 a nZ n 

then 

a " = TT~~\ (0°g(n))%(n) + . . . Qog(n))%(n)) 
where (l)j( z ) f or J = 0, . . . , /? are Gevrey-1 series with rational coefficients. 

Recall that a series X)^Lo an2; ™ ^ s Gevrey-1 (resp. arithmetic Gevrey-1) if ^2^ =0 {a n / 'n\) z n is convergent 
at z — (resp. a G-function). Lemma 12.41 and a deformation of the contour argument implies the following. 
See also [CGl Sec.7]. 

Proposition 2.5. If G(z) = ^2^ =0 a n z n is a G-function, then (a n ) is of Nilsson type, where A in Equation 
^ are the singularities ofG(z), and a, (3, f\, a ,p are determined by the local monodromy of G(z) at z = A. 

Remark 2.6. In fact, Jungen's proof combined with Andre's theorem that the series h a ,p{w) are G-functions, 
implies that the series f\, a ,p{z) of Equation Q are arithmetic Gevrey-1. 
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2.4. Proof of Theorem \5[ The next lemma is well-known. 

Lemma 2.7. If (a n ) is holonomic, the generating series G{z) = Xl^Lo ™- 2 ™ * s holonomic. 

Birkholf-Trjitzinsky, followed by Turrittin (see |BT[ ITui IPrj and |BC[ Eqn.1.3]) prove the following result 
concerning the asymptotic expansion of a holonomic sequence. 

Proposition 2.8. If G(z) — X^J^Lo ™ -2 ™ * s a holonomic function, then 

(13) a„~ «! s A-"n- Q - 1 (log(n))' 3 / A , a , /3 , s 

where A lies in a subset of the finite set of singularities of G(z), s lies in a finite set of nonpositive rational 
numbers, and a, (3 are the exponents in the local expansion of G(z) around X, and f\,a,p,s(z) « r e Gevrey-1. 

2.5. The exponents of the sequence of Theorem [TJ It remains to compute the exponents of the 
holonomic function G(z) associated to the sequence of Equation ©. One way to solve this problem is to 
convert the holonomic equation ((Sj) into a differential equation for the generating series and compute the 
exponents of the differential equation using Frobenius's method; see [Oj|In]. In addition, one needs to show 
that the corresponding constants c a ,p in (| 1 2jl are nonvanishing. An alternative way is to relate the exponents 
of the generating series G(z) of a sequence (a n ) to the asymptotic expansion of the sequence itself. 

Consider the sequence (o„) given by §8§ and its generating series G(z). Converting the recursion relation 
for (a n ) into a differential equation for G(z) we obtain that G(z) satisfies the inhomogeneous differential 
equation: 

(14) z(z 2 -7z + 2)G'(z) + {z 2 -4z-3)G(z) + z = 0, G(0) = 

If we wish, we can divide by z and differentiate once to get a linear second order differential equation for 
G(z). The singularities A of G(z) is a subset of the roots of z(z 2 — 7z + 2). I.e., we have: 

(15) Ac{0,i(7±V33)}. 
Frobenius's method gives that the exponent at A± = |(7± \/33) is given by 




do) °^ = - 1±5 2)fri 

which is non-rational. It is easy to compute that (3 — 0. It remains to argue that the so-called Stokes 
constant c Q±i/ g 7^ 0. One can do an explicit numerical computation in the spirit of |FT| Sec. 4], using Fade 
approximants and working in the so-called Borel plane. 

Alternatively, we may argue as follows. If G(z) is analytic at j(7 + \/33), then by Galois invariance and 
Equation (p~5|) . it follows that G(z) is entire. Lemma [2~9l below implies that G(z) is a polynomial. It follows 
that a n = for sufficiently large n. The recursion relation ([8]) implies that a n = for all n £ N, an obvious 
contradiction. □ 

The next lemma was communicated to us by Y. Andre, and is a useful way of deducing the existence of 
singularities of G- functions. For a detailed discussion, see also |C-L] . 

Lemma 2.9. Every entire G-function is a polynomial. 

Proof. According to Chudnovsky and Katz, a G-function is a solution of a Fuchsian differential equation, 
i.e., regular singular in C U {00}. An entire G-function does not have any monodromy at finite distance, 
hence it does not have any monodromy at infinity as well. According to a classical result of Schlesinger, 
any solution of a Fuchsian differential equation which is invariant under the global monodromy group is a 
rational function. If, moreover, it is entire, then it is a polynomial function. □ 
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3. Further discussion 

Theorem [5] may be viewed as a way of constructing holonomic G-functions from enumerative combi- 
natorics. There are two well-known sources of G-functions: from arithmetic (see Theorem [6] and also 
[An[ IBo( [DCS, Kaj), and from geometry, related to the regularity of the Gauss-Manin connection. For 
the latter, see for example example, }Br | IDe|lKa] . In all cases (combinatorics, geometry and arithmetic), the 
constructed G-functions are regular holonomic with rational exponents. 

The G-functions obtained geometry and arithmetic are closely related. The main conjecture is that 
all G-functions come from geometry. For a discussion of this topic, and for a precise formulation of the 
Bombieri-Dwork Conjecture, see the survey papers of |Bo[ IKaj and also [Tot p. 8]. Our question is motivated 
by Theorem [5] and Bombieri-Dwork Conjecture of |Toi p. 8]. 

Question 1. If (a n ) is an integer valued, exponentially bounded holonomic sequence, does it follow that it 
is a multisum sequence? 

Our next question compares the G-functions of Theorem [2] with those that come from geometry. 
Question 2. Does every G-function of Theorem [2] come from geometry? 

In [Gal] this was shown to be true when the balanced term t„ ^ is special, i.e., it is a product of binomials 
of linear forms of (n, k) (in other words, ej = +1 for all j = 1, . . . , J in Equation ©. 

Finally, let us point out that the proof of Theorem [T] and Theorem [5] is not constructive. In particular, 
it would be nice to be able to compute the singularities of the generating series of a balanced multisum 
sequence directly from the balanced termt. With this in mind, the author developed an efficient ansatz for 
the asymptotics of balanced multisum sequences; see |Galj . When r = 1 in Equation |2]) (i.e., for single- 
sums), the ansatz can be proven using the Euler-MacLaurin formula and various ideas of resurgence; see 
[Ga2] , 
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